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ABSTRACT 

This paper reports an investigation of the dynamic behavior 
of a pin-ended inelastic column subjected to a half-sine dis- 
placement pulse at one end applied with small eccentricity, 
The rates of end displacement are such that the effects of 
axial inertia are considered negligible, The column is re. 
placed by a lumped parameter mathematical model and the equa. 
tions for the model are solved by a high speed digital com- 
puter, The stress-strain relation during loading is assumed 
to be of the Ramberg-Osgood type. Unloading occurs along a 
line with slope equal to that of the tangent at the origin, 
Slendermness ratios of 50, 70, and 100 and two different values 
of eccentricity are considered in the investigation, The max. 
imum load supported by the column and the residual lateral 
deflection at tne center of the column are found as functions 
of the displacement pulse heignt and of the displacement pulse 
duration, 

It is snown tnat an inelastic column, like the elastic 
column, can support dynamic loads considerably in excess of 
the static load capacity. It is also shown tnat witnin the 
range of pulse durations studied, the smallest residual de- 


flections generally accompany the largest maximum loads, 
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In the study of structures the column ranks as one of 
the more provoking elements as judged from the amount of 
literature published on its behavior, The development of 
elastic column theory is comparatively advanced and dates 
back to 1757 when Euler published his analysis of the strength 
of columns, In more recent years inelastic column theory has 
been developed to tne point where a satisfactory explanation 
has been obtained for the phenomenon. Shanley /1/* nas been 
credited with revising and reinterpreting inelastic column 
theory from previously accepted erroneous conjectures, Util- 
izing studies of a highly idealized column he concluded that 
the tangent modulus load is the maximum load at which an in- 
itially straight, centrally loaded column will remain straight; 
loading beyond this load causes pending waich produces per- 
manent deformation, 

Tne dynamic behavior of structures has become an increas- 
ingly important engineering consideration in recent years and, 
with the aid of more sophisticated devices for the problem 
solution, many more phenomena can de included in an analysis 
than could previously be reasonably considered, The behavior 
of dynamically loaded elastic columns has received consider- 
able attention and a number of investigations nave resulted, 

Hoff and associates /2,3/ considered the case of an elas- 


hie column of uniform cross section, initially curved in the 


# Numbers in / / refer to entries in the bioliographny 





shape of a half-sine wave, subjected to constant velocity 
loading such as that encountered in compression tests in 
commercial testing macnines, It was shown that the maximum 
loads supported by the column are greater than the Euler load 
when the loading is rapid, the column is slender and the ini- 
tial deflections are small. Sevin /4/ confirmed the previ- 
ous results while including the effects of axial inertia 
which were not considered by Hoff. He inferred that, for an 
elastic column, axial inertia effects are of negligible im- 
portance with regard to the gross benavior of conventional 
structural columns, regardiess of the initial deflected shape, 
end constraint, or type of axial loading. 

Taylor /5/ investigated the case of an eccentrically 
loaded elastic column subjected to a half-sine force pulse 
and ne included the effects of axial inertia, Using a limit- 
ing value of extreme fiber strain as the failure criterion 
ne jJetermined the column fallure load quantitatively. He 
demonstrated that elastic columns will support loads much 
greater than the Euler load when subjected to rapid loading 
and that axial inertia effects become significant with this 
type of loading as the rapidity of loading 1s increased, 

Gerard and Becker /6/ studied column behavior under con- 
ditions of impact, by utilizing the properties of a propa- 
gated elastic compressive stress wave produced on failure 
of a tension specimen. They concluded that a column can 


momentarily support a dynamic compressive stress of any 





magnitude and that buckling may be confined to a length less 
than the entire length of tne column, 

The inelastic column was investigated by Chawla /7/ 
in essentially the same manner as Hoff (with an initial 
curvature or effective eccentricity) and ne found that for 
relatively slow rates of end displacezent th? maximum column 
load was less than the tangent modulus load. Brooks and 
Wilder /8/ studied the effects of dynamic loading on an 
idealized H-section inelastic column with an initial lateral 
deflection using a constant velocity end displacement. From 
their results they concluded that the static maximum load 
may oe employed as a conservative estimate of the maximum 
joad of a column regardless of the rate of end displacement, 

Investigations of the dynamic behavior of inelastic 
columns are limited in number and the problem of determin- 
ing the benavior of inelastic columns subjected to pulse 
loading appears to have been neglected. Whereas constant 
yelocity loading leads to the ultimate Tailure of the column, 
pulse loading may not nscessarily result in railure., Pulse 
loading allows the column to recover from a deflected rposi- 
tion with an amount of residual deflection while retaining 
& capacity to sustain subsequent load, The problem con- 
sidered in this investigation is that of determining the 
dynamic behavior of a pin-ended inelastic column subdjected 
to a half-sine displacement pulse at one end applied with 


small eccentricity. The problem involves the prediction 





of the maximum force developed in the column for a specified 
pulse neight and duration of loading. The problem also in- 
eluijies the prediction of the residual lateral deflection 
resulting from the pulse loading, The effects of pulse height 
and duration are investigated for two valuss of eccentricity 


ana three column slenderness ratios, 
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MODEL. 

The column considered in the investigation is initially 
straight, pin-ended, and of constant cross-section, The ef- 
fects of axial inertia and shear strains are neglected, but 
non-linear axial strain components due to bending at the 
center portion of the column length are considered, The 
effect of a non-linear stress-strain curve on Joading and 
departure from this curve on unloading is also included, 

The loading on the column is produced by axial displace. 
ment of the pin at one end, The loaded end is free to rotate 
&bout the pin, which is eccentrically placed, The other end 
is constrained by a fixed pin of equal eccentricity. The 
end displacement is imposed by a moving ram executing a half- 
sine displacement pulse. The ram can exert only compressive 
force on the column (see Fig . 1). 

The problem is developed and solved, not in terms of the 
real column, but in terms of the lumped parameter mathematical 
model shown in Fig. 1. The equations for the model are then 
solved in a high speed digital computer, The force-deflection 
characteristic of the springs is derived from the material 
stress-strain diagram, The assumed behavior 1s indicated in 
Fig. 2 where the solid curve OABD represents the variation 
of compressive stress © under montonically increasing come 


pressive strain € . If, upon reaching roint B, the strain- 
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Fig. 1 
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irg is reversed, the path will be along BC whose slope is 
that of the tangent to the solid curve at G. A second re- 
versal will cause the compressive stress to increase again 
following CB. After again reaching the curve at B, the patn 
will be along the solid curve tovard D, Additional strain 
reversals would result in additional "detours" parallel 

to BC, In the study the stress-strain relationship after 
departure from the compressive stress-strain curve is con- 
sidered Jinear into the tensile stress region. The maximum 
tensile stresses in the column model are limited to values 
near the nominal yield stress of the material. It is pre- 
sumed that any tensile stress in the column model springs 
is insufficient to cause plastic straining, 

This behavior is incorporated in the springs by taking 
the stiffness of each as K=A,E,/), where A, is half the col- 
umn cross-sectional area, E; is the tangent modulus, and X 
is the unstrained length of a spring. Again referring to 
Fig. 2, at any point on the solid curve E, is the slope of 
the curve, For points below the curve E-=E (Young's Modulus), 

This model having only two degrees of freedom was cnosen 
for the investigation in order to limit the computer solution 
time to reasonable values, but yet to retain the effects of 
Inelastic column behavior, It is felt that this model is a 
good compromise in that it aprroaches the deflected shape 
of the inelastic column and also retains the effects of 


lateral inertia and of the longitudinal strains on each side 





of the column due to both the axial load and bending, 
From the free body diagram of half the column model 
shown in Fig. 3, using d'Alembert's principle, the follow- 


ing equations may be written: 


P=FAi+h (1) 
q=m oy (2) 


P(e cos@ + + sin@)— Q(> cose ~-e sin) = 


Dividing the latter by cos@ and substituting for @ from 
equation (2) 


P(e += Tane) — im. dy (+ — etane)= 


(3) 





Fig. 3. Free Body Diagram of Half the Column Model, 








From kinematic considerations of the column model in 


a deflected position the following equations may be written; 


d= db, +2E SiN (4) 


do= dm —- 22 Sine (5) 


where <, and cf, are the axial shortening of Spring 1 and 
Spring 2 respectively and Sin is the average axial short- 


ening at the center of the model, Further 


K=dm+2:E(i-cose) +2 e sine (6) 
alone 
Yisee SS we(i=cece) (7) 


where x is the end displacement at tue loaded end of the 
model and y is the lateral deflection at midlength in the 
model. 

The following approximations are consequences of assum- 
ing that © is small and e/L is small; 


= >> etane 


as Sin@ >> e(\-—coss) 


=) 


¢ 


Cane wv sine 


fy) 


Oe 
|-cose ~ 7? 


An evaluation of the errors resulting from tnese arproxima- 





tions is given in App. VI. 
Equation (7) now becomes: 
- 
° ae ad (8) 
Now, substituting the approximations, together with equa- 


tion (&), into equations (3) through (6) and rearranging, 
results in 


aoe = Hr (ety) — (i-Fa) BE 





(9) 
Jf, = Sn a 


rb 


P| 


$=d, -t2y 


X= bn + = yo+ 48 y 


ie 


solving for dyn and substituting in the above expressions 
for of, and =A gives 


J, 


-— —&. y* e 
AY ABY + ARY 


(10) 
apex — ay” — fey —-38y (11) 


Since the model has but two degrees of freedom, only a 


limited correspondence may be established between column and 
model perameters, 


Energy considerations are used for this 
purpose, 
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To determine this correspondence both the column and the 
model are considered for the elastic case, Figure 4 indicates 
an elastic column in a deflected position during free vibra- 
tion. It has an assumed half-sine deflected shape and the 
motion 1s considered sinusoidal with time, 

The maxinum potential energy is found from the elastic 


strain energy due to bending and is 





He ay 
Potential Energy = a ( (g Z| AX 


° 


Oo 
x 
Pr 


ee cel” 
where aA Se. 


is the moment of inertia of the column cross-sectional area, 


48 the curvature of the axis of the column and [I 


I is equal to en A, where A is the column cross-sectional 


area, This can be shown to be 


2 
Potential Energy = 6 a (column) 


The maximum kinetic energy of the column, neglecting 


axial inertia effects, is 
" 2d 
Kinetic Energy =3 Kr \ yar 
0 
which can be shown to be 
2,2 |_ 
Kinetic Energy =2 HW Wr ym (column ) 


where pA 1s the mass per unit length, UW, 1s the circular 
freauency of first mode lateral vibration of the column and 
L is the column length, Setting the maximum potential energy 


equal to the maximum kinetic energy and solving for the circular 


dig: 








| 


cre 


Tx 


y = Yn Bin L 


Fig. 4. Elastic Column in Free Vibration. 





Fig. 5. Elastic Column Model in Free Vibration, 





Ol feaii- .» 


frequency of first mode lateral vibration of the column re- 


sults in 


Go eye ae f+, (column) 


Figure 5 shows the model with non-eccentric pin connec- 
tions in a deflected position, The stiffness of each spring 
is Si and the motion of the mass is again considered sinusoi- 
dal with time. The maximum potential energy as found from 


the elastic strain energy in the springs 1s 
ri 2 
Potential Energy = 4 AEC S, + S. ) 


Considering the model to be in free vibration without axial 
load the average spring deflection is zero or & =- 5 A 
Under these conditions it can be shown from equations (10) 


and (11) that 
2. 2 ce 
= _ ie é 
Sad. = 1% D/ 
and the model maximum potential energy becomes 


a 2. 
Potential Energy = rs ee A (model ) 


The maximum kinetic energy of the model neglecting axial 


inertia effects is 
2a 
Kinetic Energy = 4MWry, ym (model) 
Now, setting the maximum potential energy equal to tne max- 


imum kinetic energy and solving for the circular frequency 


a) 





of first mode lateral vibration of the model results in 


Un = / BER (model ) 


To determine the undeflected length & of the model 
springs the potential energies of the column and model are 
equated. Taking the center deflections y,, of column and 


mogei te be equal. X is found to be 


R= 26 (12) 


To determine the relation between column and model 
masses the kinetic energies of tne column and model are 
equated. Again taking the center deflections y,, to be equal, 


mis found to be 


L 
i ao (13) 


With the two previous correspondences it can be shown that 
the circular frequencies of the column and the model are the 
same, 

Now, considering the differential equation of motion 
(as previously derived for the model) for the elastic case, 
Copether with the pesults of the elastic enersy developiient, 
Euler Joad and end displacement corresrondences bdetween the 
column and model are determined, These will be used later in 
the analysis to define non-dimensional losd and disrlacement 


paramevers, If the acceleration and the eccentricity or the 





column model are both zero, equation (9) becomes 
Py =(F,-F2)e (14) 


Equation (14) applies to static buckling of an elastic column 


With zero eccentricity. Now, for the elastic case 


i= AE ma in HE, 
or A 
poh AE (d,~da)= 3 AE (S-d) 


From equations (10) and (11) 
J,- d= 5&8 y 
Thus a4 oAE 
rioPea a> Oe y 
Substituting this into equation (14), the Euler load P, of 


the column model is found to be 





where 1/Q is the slenderness ratio of the column. In terms 


of the total cross-sectional area the equation becomes 


ee 
Me Swe)" (15) 


The Euler load, in terms of the Euler displacement, is 
be — aK AE 


and, substituting for K and P, , the Euler displacement becomes 


AO 


Kem mh (16) 


15 








With these results determined from elastic behavior of 
the column and model, the equations for tne strains and for 
dynamic equilibrium of the inelastic column model are develop- 
ed in a convenient non-dimensional parametric form, From 


equation (10) the strain in Spring 1 is 


c eh - 
ye $5 


or 


é- ae | * -4%y* 4ey+4Ry | 


ios oes ee 
Defining a non-dimensional lateral deflection as V= y/e and 
making the substitution y=Ve, together with the substitution 
for the Euler displacement from equation (16), the strain may 


be expressed as 
mW. | X _ pete 2 
Geet | TEST CEIV+ EV] an 


Applying this same reasoning for the strain in Spring e re- 


sults in 


é,= aR e- o- ($) Vv - $v (18) 


Xe 


From these equations for the strains at the center for- 
tion of the column model it can be recognized that the first 
terms are the dimensionless end displacements, the second 
terms are the corrections to this end displacement due to 
the angular rotation of the rigid rods, the third terms are 
the corrections due to the change of the effective eccentri- 


city with lateral deflection, and the fourth terms result 


16 





from bending. 

The relation between the strains and stresses during 
initial loading of the column model is given by the Ramberg- 
Csgood /9/ equation: 

~ | n 
=£ 42(EY (£) 

Ce +4 (Ss E 
where € is the compressive strain, 0 is the compressive 
stress, Oy is the yield stress of the material and n is 
the shape parameter of the stress-strain curve. By defining 
a non-dimensional stress S as icon the previous equation may 


E 
be written ag 


e=S+Fbg,) S (19) 


This is the equation which is used to determine the stress 
in each spring as the model is initially loaded, For the conm- 
pressive stress-strain relation considered in this investi- 
gation the shape parameter is n=10 and the dimensionless 
yield stress is Sy=0.004, The stress-atrain curve is shown 
in Fig. 6 and is typical of 245-T aluminum alloy material, 
As each of the springs unloads, it leaves the stress- 
strain curve following a straight line of slope ¢> = leas 
As subsequent loading takes place in either spring the same 
straight line is retraced and, as the stress reacnes a level 
that puts it on the stress-strain curve again, it follows 
the curve until further unloading occurs, 


The differential equation of motion is developed in 


ee 
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nondimensional form by dividing equation (9) by the eccentri- 
city e and the square of the circular frequency of first mode 
lateral vibration Gor and substituting from equation (1) for 


FP, After rearranging terms the equation becomes 


1 de) ~ 4 (Lo) (hak )- F=f | 
ea #, (-/e) ores e+ AE "a 
A dimensionless time T is defined as T =W),t, so that 

Pa a 
Woe (at)> = (aT). Then -b¢ aes a becomes oie) Ys) 
n 
> 
By making the substitution V for y/e this becomes 2 which 


a= 
is symbolized as V. By substituting V and V and recognizing 











tne form as the previously defined non-dimensional stress, 


= 
AE 
the differential equation of motion becomes 


Vee (LR) (S.4 S,)(V+ )-(S,-S)&| (20) 


The two terms in the brackets can be recognized, respective- 
ly, as proportional to the external berding moment due to the 
force on the end of the column, and the internal bending mom- 
ent due to the spring forces, 

In order to solve this equation it 1s necessary to specify 
the forcing function, This is accomplished by specifyirig the 
displacement of the ram (as shown in Fig. 1) as a function of 
time, When the ram is in contact with the column model the 
ram displacement is the end displacement of the column, If 
the ram is not. in cortact with the end of the column, the 


force is zero (S,=-S,) and the column exhibits free vibration, 





In either case the conditions of compatibility and dynamic 
equilibrium are maintained, 

The ram displacement is a half-aine pulse of height X 
and duration Us. To enable a selection of dimensionless 
parameters having a relationship to the column the follow- 
ing dimensionless quantities are defined: 


(a) the height of the half-sine ram displacement pulse 
as a fraction of the Euler end displacement; 


ye 
Ke 


(b) the ratio of the natural period “UU, of first mode 
lateral vibration of the column to the duration Te 
of the half-sine pulse; 


@ = 

Us 
With these dimensionless quantities and the previously defined 
dimensionless time, The dimensionless ram displacement 


becomes 
=o sing OT for of T =e 
and (21) 


oso for T > =il 


Q 


To determine the maximum force and residual lateral 
deflection resulting from a spscified ram displacement it is 
necessary to compute the dynamic response of the smodel to tne 
loading. This is accomplished through a numerical integration 
process applied to the differential equation of motion and a 
numerical iteration 3cneme to satisfy tne relations among tne 


strains, stresses, end displacement and lateral deflection, 
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3. THE NUMERICAL SOLUTION 

The equations for the model are solved utilizing FOR- 
TRAN programming and a Control Data Corporation 1604 high 
speed digital computer, The investigation is accomplished 
with two basic programs, COL 1 and COL 2, The former is em- 
ployed primarily to obtain graphical output directly to study 
the dynamic response of the model to a specified loading. Pro- 
gram COL 2 is employed to determine the numerical values of 
maximum force and residual lateral deflection for a specified 
loading on the model, The complete programs in FORTRAN lan- 
guage, together with the program notation, are given in Appen- 
ax IV. 

The numerical integration procedure is the same for 
both programs, but the program, COL 2, employed to determine 
maximum values is more comprehensive, For this reason the de- 
tailed discussion is limited to this program. A discussion 
of the differences between these programs follows the detail- 
ed discussion. Figure 7 is a simplified block diagram of 
Program COL 2 and Fig. 8 is a simplified block diagram of the 
subroutine used in both programs to determine the stress from 
the strain, 

Referring to Fig. 7, the input to the program includes 
the following dimensionless parameters; 


(a) L/e (slenderness ratio) and e/e (eccentricity ratio) 
to prescribe the column model. 


(b) c<. (ram pulse height) and @G (ratio of natural period 


of first mode lateral vibration to ram pulse duration) 
to prescribe the forcing function on the model, 


al 





Fig. 7. Simplified Block Diagram of Computer Program COL 2 


Initially SWITCH=0 
Start 


Read Input 
Output iL/e ,e/0,o, 8 on,Su, AT 
Set Initial Bo] Pe (21) pSimax 1) 2. ema es 
Conditions 
Yr/Q(T4) , WORK ,ENERGY 
Xe Passes 
Prediction Loor 
Call Subroutine 
Stresa 
Compute ¢€, 
Call Subroutine 
Stress 
Compute Vv 
Compute V, Compute yr/e 
Compute P,/P, 
Compute WORK 


Input :L/e,e/e,%,G ,n,3y,4T 















VA 


Compute ENERGY 
Write Output 


Stop 


T=T + AT 
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(c) n (Ramberg-Osgood shape parameter) and Sy (yield 
stress of the material) to prescribe the compressive 
stress-strain curve, 


(a4) AT (time increment ) to prescribe the incremental 
values in the numerical integration process, 


The initial conditions are then set to start the numerical 
integration process, 

Henceforth, the dimensional quantities and the dimension- 
less Quantities are referred to synonomously except where con- 
fusion may result, 1.e., lateral deflection for dimensionless 
lateral deflection, etc, 

The numerical integration scheme used in the programs 
is a Runge-Kutta process /10/. This process predicts the 
lateral deflection V and the lateral velocity V at a time 
T+ AT from four estimates of the lateral acceleration V cal- 
culated in the time increment AT, Accuracy is of the order 
of (AT), Within this predictor loop further calculations 
are made to satisfy the relations among the strains, stresses, 
end displacement and lateral deflection, 

The ram displacement is computed from equation (21). 

The values of ram displacement and of lateral deflection are 
then used to determine the strain in Spring 1 from equation 
(17). The main program then calls Subroutine STRESS. Because 
this routine is written with dummy variable inputs to accommo- 
date the strains and stresses in both Spring 1 and Spring e2, 
the following discussion is presented based on a general 


strain € and a general dimensionless stress S, 
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Subroutine STRESS, shown in Fig. 8, is entered with the 
present value of strain € and the previous values of stress 5a 
and strain €, calculated for the preceding time T-AT. The 
solid curve in Fig. 9 represents the dimensionless compres- 
sive stress-strain curve. The dashed vertical lines indicat- 
ed by A, B, and C in the figure are representative of con- 
ditions which must be considered in arriving at a value for 
the stress S from a value of the strain €. 

Referring now to Figs. 8 and 9, if the present value of 
strain is negative this condition is represented by A, The 
Spring has unloaded from the stress-strain curve following 
a straight line with a slope equal to unity and the stress 


4s calculated from 
S=Sa+(€ -€a) (22) 


This equation 1s also used to determine the stress S if €=0, 
If € is positive there are three possibilities to be consider- 
ed, These are; 


(a) loading is occurring along the compressive stress- 
Strain Curve: 


(b) unloading is occurring on a unity slope line; 

(c) loading is occurring on a unity slope line. 
To determine wnich of these possibilities exists the incre- 
mental change in strain, AG=E€-€, , 1s calculated, If 
this quantity is negative this condition is represented by 
B, The spring is unloading on a unity slore line and the 


stress is calculated from equation (22). If the incremental 


O4 





Fig. 8. Simplified Block Diagram of Computer Program for 
Subroutine STRESS, 


Start 


IN 


S= S.+(€ = €,) 


ne 






Sq= Sat(€ - €,) 
Compute Sr 


Fig. 9. Stress-Strain Considerations in Computer Program 
for Subroutine STRESS. 
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change in strain 1s positive two possibilities exist, Either 
the stress is on the compressive stress-strain curve or the 
stress is increasing on a unity slope line. This situation 

is represented by C in Fig. 9. The stress level is deter. 
mined by calculating a value of stress from equation (22) 

and by calculating a value of stress from the Ramberg-Osgood 
equation (19) and then selecting the smaller of these two 
values as the stress existing in the spring. If the incre- 
mental change in strain is zero, the stress 5 is equal to Ss. 

To calculate the value of stress from the Ramberg-Osgood 
equation it is necessary to compute the real, positive root 
of the equation, This is accomplished by a Newton-Raphson 
iteration /10/ where the strain € is used as the initial 
estimate for the root SS». The stress S could have been de- 
termined from the strain € each time by calculating a value 
from equation (22) and a value from the Ramberg-Osgood equa- 
tion (19) and selecting the smaller value, However, this 
would involve an iteration for each stress calculation, thus 
increasing the computer time, Therefore, this method was not 
used, 

After S. is determined the value for the strain €, is 
calculated from equation (18). The main program again calls 
Subroutine STRESS to determine the value of Sz , 

At this point the program calculates the sum 5,+5,, 
which is directly proportional to the force on the end of 


the model, to determine if it is compressive (positive) or 
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tensile (negative). As previously discussed for the model 
configuration (Fig. 1), the ram cannot transmit a tensile 
force to the model. Therefore, if this sum is negative the 
program initiates an iterative interpolation scheme to deter- 
mine the end displacement of the model whicn will cause the 
sum S,;+5, to be zero, This involves recalculating the strains 
and then the stresses for each end displacement until the force 
is zero. The model is then in a state of free vibration, This 
free vibration continues until the ram again contacts the column 
or until the calculation terminates. 

Once these states of stress (end force either positive 
or zero) are determined, the lateral acceleration V is cal- 
culated from equation (20). After four passes are made to 
this point in the program to determine the values for the 
runge-Kutta predictions, as previously discussed, the lateral 
velocity V and the lateral deflection V are calculated 
for the end of the time interval AT. 

At this point the program computes the force wnichn can 


be shown to be 
; 2 
by = Hs (Le) (S82) 


and a continuous check is made to determine the maximum value 
of this quantity which is stored for print-out. This contin- 
uous check is also made to determine the maximum value of 

compressive stress in each of the springs. These values are 


stored for use in strain energy calculations to be made later, 


ef 








A continuous check is also made to determine the extreme stress 
(tensile or compressive) in each spring for print-out. The pro- 
gram then calculates the ram work done on the model for the pre- 
ceding time increment and adds this quantity to the previous 
work done by the ram, 
This whole process continues until the ram displacement 
pulse has terminated and the column is left in a state of 
free vibration, Thsn a continuous check is initiated to 
determine the time at which S,=S,=0., The lateral deflec- 
tion at this point is the residual deflection yr in the 
model, as the lateral acceleration, the longitudinal end 
force and the internal moment are all zero, and the residual 
strains determine the permanent deformation that has resulted. 
After the residual deflection has been calculated a 
switch is set to detour this portion of the program, Tne 
new path includes a continuous check to determine the time 
at which the lateral velocity is zero (point of extreme lateral 
deflection and zero kinetic energy). At this point the net 
energy input to each of the springs is calculated from the 
area under the stress-strain curve using the present values 
of stress and the previously stored maximum values of the com- 
pressive stress, This energy is a combination of the strain 
energy of plastic deformation and of the potential energy due 
to existing elastic strain in the model, This energy should 
be equal to the ram work done on the model, 


When the energy has been calculated the program then 


28 





prints all input quantities and the results. The results 


include the values of Py/Pe ,¥r/@ »Simax » 5 and the 


2Max 
time T that each of these occur, as vell as the work and 
energy. The times are included to complement the time plots 
from the other program as an aid in visualizing some of the 
dynamic behavior and to determine the problem time termina- 
tion of the program, 

The program employed to obtain graphical output directly 
is Frogram COL 1 (the FORTRAN language listing is given in 
Appendix IV). This program differs basically from Program 
COL 2 described in the preceding paragraphs in that it stores 
the variable quantities for each time T in array form instead 
of storing only the maximum values, These array quantities 
are then available for Subroutine GRAPH (not shown in the 
listing but available on the FORTRAN Library Tape at the U.S. 
Naval Postgraduate School, and also available from the Control 
Data Corroration), Subroutine GRAPH is designed to accept 
a@amaximum of 900 points for each variable to be plotted. Thus 
the program is written to calculate only 900 points while main- 
taining the necessary accuracy for the numerical integration 
process compatible with the time increment AT, 

In addition to the input data required in Program COL 2, 
axis scaling and multiple curve data, graph title data and 
curve label data are required for Frogram COL 1. In addi- 
tion to the graphical output the program prints the input 
quantities and the values of T, x/xz, €,, €, , 5,, 8, y/e 


and F/P.for every fifteenth point. As all of these quantities 
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have been stored as 9CO point arrays any combination of these 
variables may be plotted as a function of another, Subroutine 
GRAPH converts these 900 point arrays to input data on a mag- 
netic tape for a Control Data Corporation 160 computer which 
is programmed to provide output as input to a plotter, 

An example of one set of curves used in the investi- 
gation is shown in Appendix I, The example includes three 
graphs with two curves on each graph, The first graph is a 
plot of €, and €, versus T, The second graph is a plot of 
S, and S, versus T and the third graph is a plot of y/@ and 
F/P, versus T, Other examples of graphical output are shown 
in Appendices II end III, However, the graphs shown in 
Appendix II, which are for very slow loading approaching the 
static case, have been obtained with a modified version of the 
program, The modification includes calculating 1&00 points 
to maintain the accuracy for the numerical integration pro- 
cess, then re-storing every other point of the 1800 point 


arrays as 900 point arrays to accomodate Subroutine GRAPH. 
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4, RESULTS AND DISCUSSION 


a, Character of the response 

An understanding of the results is dependent upon an 
appreciation of the distinctive features of the loading and 
their effects upon the response. Loading 1s accomplished by 
a forced displacement of one end of the column in the form 
of a half-sine pulse, This pulse is characterized by the 
dimensionless parameters © and © » the quantity < is the 
ratio of the displacement pulse height to that end displace- 
ment whicn, in the absence of lateral deflection, corresponds 
to the Euler buckling load. The parameter © is the ratio of 
the first mode natural period of lateral vibration of the col- 
umn model to the pulse duration, Thus, measures the ampli- 
tude of the loading and © measures its rapidity. 

Consider tne very slow loading (O<<!) of a slender 
column to a maximum end displacement Devore the sulier criti-_ 
al value (cx>I!). At axial loads above the Euler limit 
Ene. straight form of the column is unstable. To initiate 
lateral deflection of the column model the end pins are 
eccentrically located, The effect of this eccentricity is 
to produce bending at loads pelow the stability limit of the 
column, The maximum force (P,., ) developed in the column 
becomes less than the Euler load (P, ). However, for very 
Slow loading the maximum force developed in the column model 


is the static maximum load, Increasing values of CX above 





that corresponding to the Euler critical value do not change 
the value of the maximum force P,. 

For increasing loading rates (increasing values of ® ) 
the lateral motion of the column is retarded by the inertia 
of the mass and any meaningful correspondence to static be- 
havior is lost, However, the maximum force developed in the 
column is related to the rate of loading. Also, under pulse 
loading the inelastic column recovers from a deflected posi- 
tion with an amount of residual deflection and retains a 
capacity to sustain further load, This capacity may be 
measured in terms of the amount of residual deflection yr 
which results from the loading. This residual deflection is 
equivalent to additional eccentricity for any subsequent load- 
ing. The maximum force P,, and the residual deflection yr appear 
to be the most important performance criteria fcr the dynamic 
behavior of the inelastic column subjected to a pulse loading, 
Dimensionless forms of these two quantities are utilized in 


the graphical form of the results, 


b, Range of variables 

Solutions for P,,/P, (maximum force) and for y-/e (residual 
deflection) have been obtained for a range of c from 0.01 to 
1.0 and for a range of © dependent on the column slenderness 
ratio, The results include solutions for three values of col- 
umn slenderness ratio (L/Q ) and two values of eccentricity 
ratio (e/Q ). For slenderress ratios of 70 and 50 the mini- 


mum pulse height considered waso<= 1, For a slendermess 


Re 





ratio of 100 the minimum pulse height considered was C= 2, 
For values of ©< less than these, inelastic action is insige 
nificant, The maximum values of O< considered are based on 
the maximum value of tensile stress occurring in the column 
model springs. For the investigation the stress-strain re- 
lationship after departure from the compressive stress-strain 
curve is considered linear into the tensile stress region, 
For this reason the maximum tensile stresses in the column 
are limited to values near the nominal yield stress of the 
column material, For slenderness ratios of 70 and 50 the max- 
imum pulse height considered was <= 2,5. For a slenderness 
ratio of 100 the maximum pulse height considered was m= 4, 
For values of c< greater than these the maximum value of ten- 
sile stress occurring in the column model is considerably in 


excess of the nominal yield stress for some rates of loading. 


ec, Maximum force 

Figure 10 is a plot of families of curves of maximum force 
P,/P. versus @ for a constant pulse height C<. The families 
include slenderness ratios of 100, 7C, and 50 for an eccentri- 
city ratio of C.05. Figure 11 is a plot of the same families 
of curves for an eccentricity ratio of C.0O1. The lower limit 
of © = ©.Cl is for slow loading rates approaching the static 
case, The upper Jimit of © = 1,074 approaching: the region 
where axial inertia effects become important in the elastic 
column as reported by Taylor /5/. As the inelastic column 


model does not include the effects of axial inertia © =e 
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is the upper limit considered in this study. 

Figures 10 and 11 show that as the rate of loading is 
increased the maximum force developed in the column increases, 
This is due to the delay in the development of lateral deflec- 
tion because of lateral inertia. These figures also indicate 
that as the rate of loading is decreased the maximum force de- 
veloped in the column approaches a limit asymptotically, In 
Fig. 10 this limit may be observed for each slenderness ratio 
as all of the curves converge to a single value at Q = 0.01. 
This limit is the static load capacity of the column for an 
eccentricity ratio of 0.05. In Fig. 11 this limit is not 
attained. The spread of points indicates that for this eccen- 
triciey (e/O = ©,01) there is considerable dynamic effect at 

© = ©.01. However, the curves of minimum CX for each slen- 

derness ratio at Q = ©.01 appear to approach a limit of the 
static load capacity of the column fcr an eccentricity ratio 
of C.01. <A comperison of these figures confirms that as the 
eccentricity increases the maximum force developed in the 


column decreases, 


ad, Slow loading 

It is possible to make some observations concerning 
the values obtained at the low end of the (5 range by con- 
sidering, the static case of a centrally loaded inelastic 
column, For this case the maximum load the column can sup- 
port is between the tangent modulus load F, and the reduced 


modulus load P, . In terms of the Euler load the dimensionless 
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tangent modulus load is 


Be 2 eee 
ne (23) 


where E+ is the tangent modulus or the local slope of the 
compressive stress-strain diagram, The tangent modulus de- 
fines the value of stress at which bending begins for the 
static case of a centrally loaded inelastic column, In 
terms of the model constants developed previously this 
critical stress is 


ile oa 
cr 8 (L-/e)* 


or in terms of dimensionless stress this becomes 
. i oe _— Ss 


The ratio of Young's modulus to the tangent modulus in terrs 


of the Ramberg-Csgood parameters is 


—— n (l-o7) Se) 
a 7! + “6m ) (25) 


where n = 10 (the shape parameter) and S, = C.0C4 (the dimen- 
sionless compressive yield stress) are the parameters used 

in the investigation. By choosing a value of S,, and sub- 
stituting, this for S in equation (25) the ratio E/E, is de- 
termined, By substituting the values of this ratio and Ser 
in equation (24) the slenderness ratio is determined, In 
this way a value of the ratio E7/= is determined correspond- 
ing to each of tne slenderness ratios 100, 70, and 50. Thia 


value of &7/s eauals P,/P. from equation (23). 








In terms of the EZKuler load the dimensionless reduced 


modulus ioad is 


eae (26) 


where Be is the reduced modulus, The reduced modulus de- 
fines a maximum value of stress which can be approached but 
not exceeded for the static case of a centrally loaded inel- 
astic column, In terms of the model constants develored pre- 


viously tnis critical stress is 


— * Ee 
Ocr ras 8 ey 


or in terms of dimerisionless stress this becomes 
eas il 9 as) 


The reduced modulus is derendent on the tangent modulus and 
the geometry of the column, For the column model it can 
be shown that the reduced modulus is 


x. Cee 
EK E+ E+ 


or in dimensionless form this becomes 





Ee a 2 Er/E 
BE + E/E (28) 


and this value is equal to P,/P. from equation (25). By 
cnoosing a value of Se, and substituting this for 5 in equa- 
tion (25) the ratio £E+/H is determined, The ratio =£,/E also 
determines =,/ from equation (28). By substituting the values 
of £e/t and S.. in equation (27) the slenderness ratio is deter- 


mined, In thts way a value of the ratio Ep/E is determined 
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corresponding to each of the slenderness ratios 100, 7C, and 
50. This value of E,/E equals P,/P. from equation (26). 
The values of the tangent modulus load and the reduced 


modulus load determined from tnese two methods are: 


for L/e = 100, P,/P = 1.000 and P,/P, = 1.000; 
for L/e = 70, P,/P, = 0.969 and P,/P, = 0.9773 
for L/e = 50, P,/P, = 0,649 and P,/P, = 0.691 . 


From Figs. 10 and 11 it is seen that in most cases the max- 
imum values at © = 0.01 are less than the tangent modulus 
load although tne spread of points for an eccentricity ratlo 
of 0.01 indicates that dynamic effects are more pronounced 
the smaller the eccentricity. For tnis reason some of tne 
maximum values are greater tnan the reduced modulus load, 

To observe furthertne effect of eccentricity on tne 
load capacity of the column at very slow loading rates,a 
study was made to determine tne behavior of tne column under 
a ramp displacement of the end at a loading rate corresrond- 
ing to that available to commercial testing macnines. In 
terms of tne parameters previously used to describe the pulse 


loading of the column the ramp displacement of the nead is 


X= EXO T 
where Cc = 1,0 and © = 0,01. Computer solution results in- 
cluding a graphical presentation of Strain 1 ( €,) and Strain 
2 ( €, ) versus force (P/P,) were obtained for the three column 
Slenderness ratios of 100, 70, and 50 for eccentricity ratios 


of C.05, 0.01, and 0,COl. The graphical results shown in App.II 
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exhibit a general resemblance to the experimental strain 
versus force curves of actual column test data displayed by 
Shanley /1/. From the results given in Table 1 it is seen 
tnat for eccentricity ratios of 0.05 and 0.01 the maximum 
load is less than the tangent modulus load for all three 
slenderness ratios considered, However, tne maximum load 
increases with a decrease in the eccentricity ratio, As 

the eccentricity ratio is decreased to 0,001 the maximum 
load is greater than the tangent modulus load but less tnan 
the reduced moduius load for a slenderness ratio of 50. For 
a slenderness ratio of 70 there is considerable dynamic effect 
and the maximum force exceeds the reduced modulus load by a 


Table I Computer and Theoretical Static Load Comparison 
(For Ramp Loading of Column Model) 





i Me lle - ~ 










Static 
P./P. 
Static 
Pr /P, 
for 
e/Q = 0,05 
Dynamic ae 
(Computer) for 
2/o = 0.01 
7, a a 
for 
| e/? = 0,001 
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small amount. For a slenderness ratio of 100 there is consider- 
able dynamic effect also but the solution oscillates about an 
average maximum value of P,/ P, = 1.9. 

The tabulated values of B./P. for an eccentricity ratio 
of 0.95 for a ramp displacement loading are in agreement with 
the maximum force values P,,/P, for pulse loading at @ = 0.01 
in Fig, 10. The strain ( €, and €,) versus force (P/F. ) curves 
for tne ramp displacement loading for an eccentricity ratio 
of 0.05 shown in App. II display no dynamic behavior, In Fig. 
gunk (6/2 = 0,01), where tne dynamic effects are more pronounced 
at & = 0.01 for pulse loading,the forces on the O<= 1 curves 
for slenderness ratios of 70 and 50 agree with these tabulated 
results, For a slenderness ratio of 100 tne smallest value of 
oc considered for pulse loading was 2 and this curve appears 
to approach a value P,/P. = 0.983 asymptotically. The strain 
versus force curves for the ramp displacement loading for an 
eccentricity ratio of 0.01 shown in App. II display an insignif- 
icant amount of dynamic behavior which does not affect the 
values of maximum force P./P. ,. Therefore, these values of 
p./P. may be considered the static load capacity of the inel- 
astic column for eccentricity ratios of 0.05 and 0,01. 

Appendix III contains a compilation of the dynamic re- 
sponse curves stress (S,; and 5, ) versus time (T) for a range 
of @ from C.1 to 1.0 and for an eccentricity ratio e/Q =O god 
This graphical compilation includes tne solutions for the min- 


imum and maximum values of pulse helght considered in tne 


4] 





study for all three slenderness ratios, Included with the 
compilation is a discussion of the more prominent features 


of the curves, 


e, Residual deflection 

Figures 12, 13, and 14 are familes of curves showing 
residual deflection yr/Q versus Q for a constant pulse 
height oc and for e/Q = 0,01 and 0,65. Figure 12 is the 
family for L/e = 100, Fig. 13 is the family for L/e = TO: 
and Fig, 14 is the family for L/Q = 50, The points shown 
on tnese curves are from the identical computer solution 
results as the points shown in Figs. 10 and 1l, 

From Fig. 12 for L/Q = 100 it is shown that in general 
the values of residual deflection increase with tne rapidity 
of loading to a maximum value then decrease abruptly. This 
abrupt decrease is because of the delay in the development of 
lateral deflection due to lateral inertia, From Fig. 13 
for 1/2 = 70 and Fig. 14 for L/e = 50 it is shown that in gen- 
eral tne values of residual deflection are essentially cone 
stant as the rapidity of loading increases until an abrupt 
decrease occurs due to the effect of lateral inertia, The 
exception for ©“ = 1 and L/e = 70 is due to the small amount 
of inelastic action that takes place for tnis pulse height 
at all values of So » The amount of residual deflection may 
be expressed in terms of the amount of inelastic action tnat 
occurs, From the graphical compilation of stress versus time 


curves in App. III and tne data compiled in App. V the amount 
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of inelastic action may be determined, For a slenderness 
ratio of 100, inelastic action in the column is limited to 
maximum compressive stress levels below the value of tne 
yield stress for all solutions considered in the investiga- 
tion, For a slenderness ratio of 50, inelastic action, in 
general, occurs at much higher maximum compressive stress 
levels-near or above the yield stress, For a slenderness 
ratio of 70, inelastic action occurs at maximum compressive 
stress levels intermediate to those for slenderness ratios 
of 100. and 50, 

A comparison of Fig, 12, 13, and 14 indicates that as 
tne slenderness ratio decreases the effect of eccentricity 
on the values of residual deflection is less pronounced, 
However, decreasing eccentricity causes a shift of tne re- 
sidual deflection curves to longer pulse durations (smaller 
values of © ). This is because dynamic behavior is more 
pronounced the smaller the eccentricity and lateral inertia 
effects become more important at slower loading rates, It 
may also be observed from these residual deflection curves 
that the maximum values of residual deflection do not necessa- 
rily occur at tne slowest loading rates, The minimum values 
of residual deflection generally occur at the fastest loading 
rates, The only exceptions are for the columns of largest 
slenderness ratio (L/@ = 100) and largest eccentricity 
(2e/Q = 0,05). However, the residual deflections appear to 
follow the same trend but are limited by the maximum value of 
the parameter © = 10(shortest pulse duration considered in the 


investigation). 
46 
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533 CONCLUSIONS 


From the results discussed in the preceding section, 


the following conclusions may be drawn for the rates of 


loading considered in the investigation: 


(a) The inelastic column model, together with a high 


(b) 


(c) 


(e) 


(f) 


speed digitel computer and the peripheral plotting 
equipment, provides a sound and reasonable method 
for the study of the dynamic behavior of inelastic 
columns, 

For very slow loading (long pulse duration) the model 
exhibits behavior in agreement with Shanley's postu- 
lates. 

As eccentricity is reduced, the lateral response is 
slowed, Accordingly, columns with very small eccen- 
tricity show significant dynamic effects at quite 
Slow loading rates, 

Within the range of pulse durations studied, the 
maximum load increases monotonically with decreas- 
ing pulse duration, 

The ratio of maximum dynamic load to static load 
capacity decreases as the slenderness ratio is 
reduced, 

Within the range of pulse durations studied, the 
smallest residual deflections generally accompany 


the shortest pulse durations, 
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APPENDIX I 
Example of Graphical Output from Program COL 1 
Graphic?) Notation 
The computed points on the graphs are the values of the 
dimensionless variables developed in the basic equations, 


The graphical notation is as follows; 


Symbol] on Graph Corresponding Symbols in 
the basic equations 

SLENDERNESS L/e2 

ECC e/? 

ALPHA Cox 

BETA & 

STRAIN 1 € 

STRAIN 2 é, 

STRESS ] 5 

STRESS 2 S. 

DEFLECTION y/e 

FORCE P/P 

TIME Jk 


Each of the curves on the graph is labeled to corresrond 

to the terms in the graph title, The x-axis is the axis 

of abscissas and the y-axis is the axis of ordinates, The 
numbers specifying the coordinates are integer numbers locat- 
ed to the right of the coordinate on the x-axis and above the 
coordinate on the y-axis, The scale quantity is the value of 
the first marked coordinate from the origin in either the 


X-direction or the y-direction expressed in powers of ten 


I~] 





to locate the decimal point. The symbol E in the scale value 
indicates that the number following is the exponential power 
of ten to be applied to the basic scale value, 1i.e., 1.OOEK+t0O1l 
is: 1 eon 110" 

Discussion of the curves 

On pages I-5 to I-7 following this discussion are shown 
curves of €,, €, , 5, , 8, , y/@ and P/F versus T. These 
curves are the solutions for the inelastic column model dy- 
namic response to a ram pulse loading with the following par- 
ameters: 

(a) slenderness ratio, L/@ = 100 

(b) eccentricity ratio, e/@ = 0,05 

(c) ram pulse height,<= 4.0 

(d) ram pulse duration in terms of Qe = 0.2, 

Te = eTT/© = 107T 
The discussion includes some of the more prominent features 
of the dynamic behavior of the column model under ram pulse 
loading, 

A comparison of the strain versus time curves with the 
stress versus time curves indicates tne amount of inelastic 
action that takes place in Spring 1 and Spring 2, In Spring 1 
the compressive stress level increases above that in Spring e 
as the amount of bending increases, The internal bending 
moment in the column model is directly proportional to 
S; - Sz. As the amount of bending increases Spring @ de- 


parts from tne compressive stress-strain curve along a straight 








line with unit slope and continues into the tensile stress 
region, The values of stress in Spring 2 remain on this 
straight line for the remainder of the solution. The con- 
pressive stress level in Spring 1 continues to increase to a 
point where it also departs from the compressive stress-strain 
curve along a straight line as the internal bending moment in- 
creases and the force decreases to zero, This force is dir- 
ectly proportional to 5 + 5... When the force becomes zero 
the end of the column has departed from the ram which exhibits 
a half-sine pulse displacement, The column model is then in 
a state of free vibration until the end again makes contact 
with the ram, The compressive stress level in Spring 1 again 
returns to values on the compressive stress-strain curve, 
This alternate departure and contact between the column and 
the ram continues until] the ram displacement pulse nears its 
termination, The column then contirues in a state of free 
vibration where S, = = &:., The amount of residual deflec- 
tion may then be determined from the mean value of the deflec- 
tion or from the deflection at the time when S. = &S- = 0, At 
this point the internal forces and the internal bending moment 
in the column model are zero, 

The effect of eccentricity on the dynamic behavior of 
the column model with all other parameters remaining the 
same may be observed by comparing the stress versus time 
curves in this appendix with the stress versus time curves 


on page III-13 where the eccentricity ratio is 0.01. It may 


I-35 


(i 








be noted “hat the sssential differences in these curves are 
that the maximum compressive strsss levels attained in Spring 1 
and Spring 2 are greater for the smaller eccentricity. Also, 
contact between the ram and the column terminates sooner for 


the smaller eccentricity. 
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APPENDIX If 


Theoretical and Computer Solution Graphical Comparisons of 


Maximum Loads for the Static Case 


The notation for the graphs displayed in this appendix 
is explained in Arp. I. 

The graphs on pages II-3 to II-ll folowing this dis- 
cussion are the curves ~, and .,versus P/F for the three 
slenderness ratios considered in the study. As previously 
mentioned, the loading for the study in this appendix is a 
ramp displacement of the end of the column model at a dis- 
placement rate comparable to that in commercial testing ma- 
chines, The study includes eccentricity ratios of 0.05, 
O,.01, and 0,001 for each of the slenderness ratios, The 
tabulated results of this study are shown in Table I in 
Section 4, 

In addition to the previous comments made concerning 
the results, it may also be seen from the graphs that for 
all slenderness and eccentricity ratios considered an amount 
of bending takes place before strain reversal occurs in 
Spring 2 and before the maximum load is reached, This was 
essentially Shaniey’s point in proving the fallacy of the 
reduced-emodulus theory for column action, As stated by 
Shanley /1/, 

oooocootvne reduced-modulus theory is nat correct 


for predicting the maximum load up to which a pere 
fect column will remain straight. This is because 


{It-l 





it is possible for the column to bend simultaneously 
with increasing axial load, Under such conditions 
it is possible to have bending without introducing 
any strain reversal, upon which the reduced=modulus 
theory depends, 
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APPENDIX III 
Column Model Dynamic Response Curves 

The notation for the graphs displayed in this appendix 
is explained in App. I. 

The graphical compilation on pages III-3 to III-26 
following this discussion givestne curves 5S, and Sz versus 
Jk for @ = O.,1, C.2, 0.5 and 1.0 and for an eccentricity ratio 
e/@ = O.Cl. Included in the compilation are the solutions 
for the minimum and maximum values of ram pulse neight CX con- 
sidered in the investigation for all three slendermess ratios, 
The graphical solutions are grouped by slenderness ratio with 
the solution for the maxinum value of <& followed by the solu- 
tion for tne minimum value of &% for a single i) value. The 
order of the & Values is C.1, 0.2, OuS @nd 1.0, that is. trop 
the longest pulse duration of loading to the shortest pulse 
duration of loading consijered in the compilation, 

A discussion of some of the prominent features of tne 
dynamic behavior of the column model under ram pulse loading 
is included in Apr. I where the curves of €,, €,, S,, Sz, y/Q 
and P/P, versus T are displayed for one computer solution, 

The following is a discussion of some of the more prominent 
features of tne dynamic benavior of tne column model in terms 
of tne stress curve compilation in this appendix. 

In general the amount of inelastic action in the column 
model is inversely related to the slenderness ratio, Fora 


Slenderness ratio of 100 inelastic action is limited to com. 


III-1 





pressive stress levels below the value of the yield stress. 
For a slendermess ratio of 50 the maximum conpressive stress 
levels are near or above the yisld stress. For a slenderness 
ratio of 7O inelastic action occurs at maximum compressive 
stress levels intermediate to tnose for slenderness ratios 

of 50 and 100. 

It may also be noted that, in general, contact is lost 
between the ram and the end of tne column (when S, = - 5,) 
before the loading pulse nears its termination for ©) values 
of ©.e and 0.5 and slenderness ratiosof 70 and 50. For a 
glenderness ratio of 100 this occurs at values of e = O75 
G22%and O55. 

It may be seen that the value of maximum tensile stress 
in Spring 2 exceeded the value of the nominal yield stress 
(Sy = C,004) for the e values of ©.2 and 0.5 and &< #2 2,5 


for a slenderness ratio of 50, 
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APPENDIX IV 


Comruter Programs and Notation 


- —~— 


ee oe deta 

The FORTRAN language computer programs, togetner with 
tue notation, and the progran operating instructions for 
bota Prozram COL 1 and Program COL e@ are included in tne 
following pages. The discussion and a listing of Prosran 
COL ] are included on pages IV -2@to IV -10. Tne discuss- 
ion and a listing of Program COL 2 are incladed on pares 
TV -ll to IV -20, The FORTRAN language format for each of 
the crograms is in standard FORTRAN language to oe used with 
a standard FORTRAN compiling routine, The program largusage 
is compatible with any of the large machine FORTRAN comriling 
systens including tne I3M 7090 and the CDC 1504 systens, Each 
line in the listing is tne information contained on one &0 col- 
umn IBM card where columns 1-5 contain any statement numboder, 
column 6 is used to indicate a continuation card, columns 
7-72 contain tne FORTRAN statement, and columns 73-&0 are 
used for card identification only. Altnousghn the notation is 
basically the same, the subscript notation is different for 
each program, Program COL 1 stores a number of varlabis 
quantities in array form for rrint-out and to accommodate 
Subroutine GRAPH, whereas, Program CCL @ stores only one 


value of tne variable quantities for print-out, 


1V-1 





Pesivevat ion for Program COL 1 

In tne notation for Program COL 1 tne subscript J re- 
fers to quantities within the Runge-Kutta prediction loop 
which are calculated for intermediate times within tne tine 
increment AT. The subscript I refers to tne quantities 
corresponding to each time T(I). The notation for Program 
COL 1 is as follows: 


Computer program notation for Program COL 1 


—-— =o 


Computer Corresponding Computer Corresponding 
Progran terms in table Program terms in tables 
_fern _of symbols _term _of symbols __ 
A n SIGY Sy 

ALPHA ot SIG1(J) SF 

BETA Q SIG2(J) sv 

DDV(J) V SLNESS L/e 
DEFL(T) V STRN1(I) S 

DELT AT STRN2 (1) ee. 
DISPL(T) Kee STRS1(T) S\ 

ECC 6/2 ST RS2(I) S, 
EPSI1(J) c, ee, I 

EPSI2(J) €. VEL V 

FORCE(I) P/P, X( J) X/Xe 


ec, Operating instructions for Program COL ] 
The output from Program COL 1 includes tne print-out 
and the array quantities necessary for Subroutine GRAPH, 


The listing for Prozram CCL 1 does not include Subroutine 


TV-2 


re 





GRAPH, However, a listing of Subroutine GRAPH and instruc- 
tions for its use may be obtained from tne Control Data Cor- 
poration. To use tne program to obtain tne print-out quanti- 
ties without tne graphical output remove tne cards CALL GRAPH. 
To operate Program COL 1 without Subroutine GRAPH use the cards 
including Subroutine STRESS in the order shown in the listing 
followed by one blank card, then the input data card, 

The read-write statements in the listing are written 
for magnetic tape input and outnut as follows: 

(1) Input (READ INPUT TAPE 2) 

The input to the program is contained on one card (7#10,.3) 


and includes the following seven quantities in order; 


Input quantity Format 
SLNESS (L/@ ) E10.3 
ALPHA (CX ) E10.3 
ECC (e/Q) HIG 
BETA (“© ) OO Pas, 
A (n) BOR 
SIGY (Sy) Z10.3 
Dane { AT) E10, 3 


If Subroutine GRAPH is used axis scaling and multiple 
curve data, graph title data and curve label data in 


addition to the above are also needed in tne input. 


(2) Output (WRITE OUTPUT TAPE 4) 
The outout includes tne input quantities and the vari- 


able quantities for every fifteenth noint as follow; 
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Outnut Quantity Format 


SLNESS (1/Q ) F5,0 
ALPHA ( co ) P52 
ECC (e/Q ) F5.> 
BETA (@ ) F4,2 
A (n) F5.1 
SIGY (Sy) foes 
DELT (AT) Shs ee 
T Lote 
DeSrin ( x7 x) £14 ,4 
Sale G)) 214.4 
STRN2 (€.) 214.4 
STRS1 (5, ) E144 
STRS2 (S2) 214.4 
DEFL (y/Q ) 214.4 
FORCE (P/P- ) alee 





PROGRAM COL 1 

C PROGRAM TC DETERMINE INELASTIC CCLUMN RESPONSE TO HALF-SINE 

C RAM DISPLACEMENT PULSE APPLIEC AT END WITH SMALL ECCENTRICITY 

C TO INITIALLY STRAIGHT CCLUMN. WHEN RAM CAUSES TENSION CCLUMN 

C EXHIBITS FREE VIBRATION. 
ODIMENSION T(900),CISPL(9CC),CEFL(900) »STRNI(9CO) »STRN2(9CO), 
ISTRS1(900),STRS2( 9OC) »X€4U) ,OOV(Y) ,EPSIV(U) ,EPSI2Z(UY,SIGI(4)» 
2S1G2(4%) »AK (4,2) 2C (4), FORCE(900) 
COMMCN A,SIGY 


READ INPUT TAPE 2¢5sSLNESS ,ALPHA,ECC, BETA, A,SIGY,DELT 


KA 


FORMAT(7E10. 3) 
TCT) =C.0 
DEFL(1)=0.0 
VEL=0.C€ 
STRN1(0)=0.0 
STRN2(0)=0.0 
STRST(C}=0.0 
STRS2(C}=0.C€ 


tlic 
KTENS=} 
KONTACT=} 
M=0 


IV-5 





I=] 
bO CO 70 J=1,h 
TC=T(1)4C(J) 2 DELT 
VC=DEFL( 1 )4+C( J) #AK(J-1,1) 
IF (0. 5*#BETA#TC—3. 141593) 50,5555 
50 X(J)SALPHA*SINF(O .5*#BETA*TC) 
GO TO 60 
55 X(J)=C.0 
6QOEPSII(J)=12.1761/( SLNESS##2 )0(X(J)-0.58( ECC #82) e VC 882+ 
LECC *¥VC#(1.0-ECC)} | 
IF (J—1)62462, 63 
62 CALL STRESS(EPSI1 (J) sSTRNICI-~1) ¢STRSICI-1V) ¢eSIGH(Y) ) 
STRNI(I)=EPSI1(1) 
STRSICI)=SIGI(V) 
GO TO 64 | 
63 CALL STRESS(EPSIV (J) ¢STRNICT Ie STRSIC IT) eSIGI(N)) 
6UOEPSI2Z(J)=12.17617( SLNESS#¥2 )e(X(JI-O.5#(ECC #82) eVCHe2— 
1ECC #*VC#(1.0+ECC)) 
IF ( J- 1965565, 66 
65 CALL STRESS(EPSI2(J5)eSTRN2ZCI-1) »STRS2(1-1),S1IG2¢J) ) 
STRN2(1)=EPSI2(1) 
STRS2(1)=SIG2(1) 
GO TO 661 


66 CALL STRESSCEPSI2Z (3) ,STRA2Z(T) eSTRS2C1),S1G205)) 


661 I[F(KTENS—1)67,665,665 


IV-6 





665 


67 


€75 


68 


685 


686 


687 


688 


689 


IF(CKONTACT—1)67,694,654 

IF (M-1)675,68,685 

X1=X( J) 

PI=SICI(J)+SIG2(5) 

X(J)=CISPL(I-1) 

M=M+] | 

GO TG 60 : 
X2=X(J) 

P2=SIGI(J)+SIG2(J) 

X( J) =X2—-( P2/(P2-P'1)) #(X2—X1) 

V=M+ ] 

GO TO 6C 

P=SIGI(J)+SIG2(J) 

IF(P)686,688, 685 
IF(10.0##(-8.0)-ABSF(P1) 687,693,693 
X1=X(J) | 

P1=P 

X(J)=xX2—(P2/(P2-P 11) 8 (X2-X1) 

GO TO 60 

IF (10.0 ##(-8.0)-ABSF(P) }689,693,693 


X2=X(J) 


- P2=P 


693 


X(J)=X2-(P2/( P2-P 1) ) 4 (X2-X1) 
GO TO €&C 
Ma 
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KTENS=1 
GO TO 698 
694 IF(SIGI(J)+SIG2(5 ))695,698,698 
695 IF(C.S#PETA#TC—3. 141593)6955, 696,696 
6955 KTENS=0 : 
GO TO 675 
696 KONTACT=C 
GO TO 675 
6980D0DV(J)=C.0410639e( SLNESS##2 )#((SIGI (J) 4S1G2 (5) # (VC+1.0)- 
T(S1G1(J)-SIG2(J)) SECC) 
AK (J, I) =DELT#(VEL+C (J) eDELT#COV(J)) 
TO AK(J,2)=CELT#O0VE J) 
DISPL(I)=X(1) 
FORCE (1 )=.0410639 #( SLNESS##2 y eC STRS ICT )#STRS2ET)) 
VEL=VEL + (AK(1,2)+2.08AK(2,2)4+2.08AK(3,2) +AK(4U 42) )/6.0 
OEFL(C I41)=DEFLUT) +(AK(1,1)42.08AK(2,1)42.08AK (3,1) ¢AK(4U,1))/6.0 
IF (1-900) 75,80, 80 
75 TCI4+))=TCI+D0ELT 
l=I +1 
GO TO &C 
8C WRITE CUTPUT TAPE 4482 
82 FORMAT(1H1) 
WRITE OUTPUT TAPE Wy GS eSLNESS ALPHA, ECC BETA, AsSIGY CELT 
8SOFORMAT( I4H OUTPUT OATA //3X 12Z2HSLENDERNESS#F5.0,2X6HALPHA=F5.2, 


12X4UHECC=E5.6 39 2XSHBETAH=FR a 29 2X ZHAZ=FE 5.19 2XSHS IG YER Ge 39 2XSHODELT HFS. 3 


= 


2 //) 
Iv-8 





WRITE OLTPUT TAPE 4,90 
SCCFORMAT(UXUHTIME 4X1 2HDISPLACEMENT p 4X1 OHDEFLECT ION, SX8HSTRAIN 1, 

T6XBHSTRESS 1,6X8HSTRAIN 2,6X8HSTRESS 2, OXSHFORCE// 1} 

DO 92 1=1,900 
92 DEFL(1)=ECC#DEFL(T) 

DO 95 I=1,9C0,I5 
QSOWRITE CUTPUT TAPE 4,98,T(1) ,OISPL(I),OEFLOI) pSTRNIC I) sSTRSIC I)» 

1STRN2Z(1),STRS2(1) pFCRCECT) 
S& FORMAT(32XF 6.2, 7ETU.4) 

CALL GRAPH(900,T, STRN1,8) 

CALL GRAPH(900, Ty STRN28) 

CALL GRAPH(SCOyT, STRS198) 

CALL GRAPH(900,T, STRS2e8) 

CALL GRAPH(900,T, OEFL,»8) 

CALL Kebenisbo.t.ebaceee 

STOP 

END 

SUBROLTINE STRESS (EPSIyEPSIA,SIGAySIG} 

COMMON A,SIGY 

[FC EPSI)20, 20521 

20 SIG=SIGA+(EPSI-EPSIA) 

GO TO 32 
21 IFLEPSI-EPSIA) 22222422 


22 SIG=SIGA+(EPSI-EPSIA) 
IV-9 





23 


24 


29 


30. 


31 


32 


35 


GO TO 32 

SIGB=SIGA+(EPSI-EPSIA) 

B=3.0/(7.0#(SIGY##(A-1.0))) 

SIGI=EPSI | 

SIG2=SIGI-(SIG14B #(SIGI##A)-EPSI)/(ABe(SIGI##(A-1.0)) 41.0) 
IFC ABSFE(SIG2#( 10. 088 (~6.0)))-ABSF(S1G2-S1G1))294 30530 
SIGI=S1G2_ 

GO TO 24 

IF (SIG2-S1G8) 31.3 1,32 

SIG=S 1G2 

GO TO 33 

SIG=SIGB 

RETURN 

END 


END 
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ad. Notation for Program COL 2 

In the notation for Program COL 2 the subscript J re- 
fers to quantities within the Runge-Kutta prediction loop 
which are calculated for intermediate times within the time 
increment AT. The variable quantities ending in B corre- 
spond to a time T - AT and the variable quantities ending 
in A correspond to a time T, The notation for Program COL 2 


is as follows; 





Computer program notation for Program COL 2 


Computer Corresponding Computer Corresponding 

Program terms in table Program terms in table 
Term of symbols _ Term of symbols 

A n SIG1(J) S, 

ALPHA oS SIG2(J) Sh 

BSTA © SLNESS L/Q 
DDV(J) V STRNLA or B & 

DEFLA or B V STRN2A or B es 

DEFLP yr /Q STRS1A or B sy 

DELT AT STRS2A or B S> 

DISPLA or B x/X¢ STRSIM pppoe Me deer 
Boo °/R srason Max Ss, (tengiue 
EPSI1(J) E, T 38 

EPSI2(J) eS VELA or B V 

PMAX Pe X(J) ye 

SIGY Sy 
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e, Operating instructions for Program COL e 

To operate Program COL 2 use the cards including Sub- 
routine STRESS in the order shown in the listing followed by 
one blank card, then the input data card. 

The read-write statements in tne listing are written for 
magnetic tape input and output as follows: 

(1) Input (READ INPUT TAPE 2) 

The input to the program is contained on one card (7£10.3) 


and includes the following seven quantities in order; 


Input Quantity rornat 
SLNESS(L/@ ) £10.3 
ALPHA (Cc ) On CO) 
ECC (e/Q) BuO 45 
BETA ( ® ) Ges, 
A (n) £10.3 
SIGY (Sy) E10 .3 
DELT (AT) NH Ore, 


(2) Output (WRITE OUTPUT TAPE 4) 
The output includes the input qQuatities and one value 
of eacn of the variable quantities including the time 


T of occurrence of tne value as follows: 


Output quantity Format 
SLNESS (L/Q) EO.e 
ALPHA (cx ) E14.4 
ECC (e/Q ) Ings Wes 
BETA (Q ) F4,2 
2. Bog. 


IVA 2 





Output guantity 
SIGY (Sy) 
DELT (AT) 


PMAX (P,,/P, ) 
Tl 

STRSIM (S, max) 
T2 


STRS2M (So max) 
tT) 


— ( yr/@) 


IV-15 


Format 
nO so 
ee, 


E14 .4 
hone 


E14 .4 
F6,2 


E14 .4 
BGe 


E14 .4 
Fe.e 





PROGRAM COL 2 
PROGRAM TO DETERMINE RESIDUAL DEFLECTION AND MAXIMUM FORCE 
IN INELASTIC COLUMN RESULTING FRCM HALF-SINE RAM DISPLACEMENT 
PULSE APPLIED AT END WITH SHALL ECCENTRICITY TO INITIALLY 
STRAIGHT COLUMN. WHEN RAM CAUSES TENSION COLUMN EXHIBITS FREE 
VIBRATION. 

ODIMENSION X(%) pODVI4) EPS ICK) EPSI2E4) pSIGI (4) SIGZC4) pAK UU 2), 

Ich) 

COMMON A,SIGY 

READ INPUT TAPE 2+5ySLNESSsALPHAyECCy BETA, AsSICY pDELT 

5 FORMAT(7E10.3) 
10 T=0.0 

DISPLA=0.0 

VELA=0.0 

DEFLA=0.6 

STRNIE=0.0 

STRN2E=0.0 

STRSIP=C.C 

STRS228=0.0 

PMAX=C.0 

STRSI¥=C.C 

STRS2¥=C.0 

STRS2E=0.0 

WORK=C.0 


Clljy=c.€ 
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C(2)=C.5 

C(3)20.5 

C(U)=1.0 

KENERG=1 

KTENS=1 

KONTACT=1 

M=0 
4O DO 70 J=1,94 

TC=T#C(J) #DELT 

VC=DEFLA+C( J) #AK( J-15 1) 

IF (0.5*BETA®TC—3. 141593) 50,55955 
50 X(J)=ALPHA#SINF(0.5#BETA®TC) 

GO TO 60 
55 x(J)=C.C 
6OOEPSII(J)=12.1761/( SLNESS##2 )#(X(J)-0.50(EGCHO2 8 VCOe2+ 

JECC#VC#(1.0-ECC)) 

LE (0-162 962563 
62 CALL STRESS(EPSII (J) »STRNIBySTRSIBsSIGI(J)) 

STRNIA=EPSI1(1) 

STRSIA=SIGIO1) 

GO To 64 
63 CALL STRESS(EPSII(J)¢STRNIAySTRSIAsSIGI(U)) 
GUOEPSI2Z(J)=12.1761/( SLNESS##2 )#(X(J)-0.5e(ECC#02)#VCee2- 
_ VECC#VC#(1.04ECC)) | 


IF (J-1)65,65, 66 
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65 


66 
661 
665 

67 


675 


68 


685 


586 


687 


CALL STRESSCEPSI2 (J) ¢STRN2B,STRS2B,SIG2(J)) 
STRN2ZA=EPSI2(1) 

STRS2A=SIG2(1) 

GO TO 66) 

CALL STRESSCEPSI2 (J) sSTRNZA,STRS2A,SIG2(J)) 
IF(KTENS-1)67,665,665 

IF (KON TACT—1) 67 9694 ,694 

IF (M-—1)675, 68,685 

X1=X(J) 

PI=SIGI(J)#S1G2UJ) 

X(J)=CISPLB 

M=M+] 

GO TO 6C¢ 

X2=X(J) 

P2=SITEI(J)I+SIG2(J) 
X(J)=X2-(P2/(P2-P1)) #(X2-X1) 

M=M+i | 

GO TO 60 

P=SIGI(J)+S1G2{ J) 

IF(P}€8@6, 688,688 

TFC 1TO.Cw# (-8.CI-ABSF(P)) 68756935693 
X1=X( J) 

Pl=P 

X( J) =xX2—-CP2/UP2-P 1) eCK2-X1} 


GO YO 60 
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688 IF(10.0##(-8.0)-ABSF(P))&89,693,693 
689 X2=X(J) 
P2=p 
X(J)=X2-(P2/( P2=-P1)) # (X2=-X1) 
GO TO éC 
693 M=0 
KTENS=1 
GO TO 698 
694 TF(SIGI(J)+S1G2(J))695,698,698 
695 IF (0.S#BETA#TC—3. 141593) 695556965696 
6955 KTENS=0 
GO TO 675 
696 KONTACT=C 
GO TO 675 
69B0DDV(J)=0.04 106398 ( SLNESS##2 )#((SIGI(J)4S1G2() )#(VC41.C)- 
1(SIGI(JI-SIG2(J)) JECC) 
AK (J, 1)=CELT#(VEL A+C (J) ®DELT#O00V(J)) 
70 AK(J,Z)=DELT*®OOVE J) 
1F(0.5#BETA#T—3.141593)7C2,704, 704 
7O20WORK=hORK+0.5#((S TRS TA+STRS2A) 4+(STRSIB*¢STRS2B) )® 
1(ALPHA®(SINF(0.5#BETA#T)-SINF(0.5#BETA®(T-DELT)))) 
704 LF(PMAX=C.0410639"( SLNESS##2 )e(STRSTA+STRS2ZA) IT) 972972 
71 PMAX=C.CH10639#( SLNESS##2 )#(STRSTA+STRS2A) 
T1=T 


72 YECSTRSIM-STRS IA) 73, 7k, Ph 
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(3 


Ty 


75 


151 


it a2 


fo5 


T56 


ot 


758 


759 


T6 


U7? 


78 


Ae) 


&0 


STRSIM=STRSIA 

STRSTE=STRSIM 

T2=T 

IF (ABSF(STRS2M)-ABSF(STRS2A) 7597519751 
STRS2N=STRS2A 

T3=T 

[LF (STRS2E-STRS2A) 75247559755 

STRS2E=STRS2A 

IF (KENERG-1)75647 6,76 

IF(VELB)758,759,757 

IF (VELA)759,759,8.1 

IF (VELA)81,759, 759 

B=3.0/(7.0#(SIGY*#(A-1.0))) 

ENERGI= (A®B/(A+1.0)) #STRSTE#*(At1.0) #0.5#STRS1B 002 
ENERG2= (A®B/(A41.0)) #STRS2E#e(A41.0) #0. S#STRS2B 022 
ENERGT=0.0821278#( SLNESS##2 )#(ENERG14+ENERG2) 

GO TO 84 

IF (KONTACT-1)77 781,81 

BENDM1=STRS1IB-STRS28 

BENDM2=STRSITA-STRS2A 

[ECSTRS1B)78,78579 

IF (STRS14)81,80,80 

[F(STRS1A)80, 80,81 

DEFLP=(DEFLA —BENCM2/ (BENDM2-BENOM1)©(DEFLA-DEFLB) )*ECC 


THT 


pai gems 





81 


KENERG=0 

T=T+DELT 

DISPLB=DISPLA 

DISPLA=X(1) 

VELB=VELA 
VELA=VELB+(AK(142)42.-08AK(2,2)42.08AK (3,2) 4AK(4,2) 1/620 
DEFLB=DEFLA 

DEFLA = DEFLB 4(AK(1,1)42.00AK( 291420 OMAK(3y 1) 4AK (491) 9/600 
STRNIB=STRNIA 

STRN2ZE=STRN2A 

STRSTB=STRSIA 

STRS2B=STRS2A 


GO TO &O 


84 WRITE OUTPUT TAPE 4,85,SLNESS, ECC »BETAsA,SIGY,OELT 


B5SOFORMAT(///14H OUTPUT DATA //12HSLENDERNESS#F5.0,2X6H ECC=F5.2, 


12X5SHBETA=FU 62, 2X2HA=F 50 1 p 2XSHSIGY=F6. 35 2XSHDELT 2F5.3//) 


WRITE OUTPUT TAPE &,86 


B6OFORMAT(6X5HAL PHA, 7X1 1IHOEFORMATION, 3X9HMAX FORCE, 5XT2HMAX STRESS 1, 


87 


88 


885 


887 


T2X12HMAX STRESS 2/7/) 


WRITE OUTPUT TAPE &,88,ALPHA, DEFLP,PMAX,STRSIM,STRS2M 

FORMAT(S5E 14.4) 

WRITE GUTPUT TAPE &,885,TH,T1,T2,73 , 

FORMAT (20X IH( F602 91H) pOX1THI( Fb 29 1H) p6X1TH( F602, 1H)» OXIK(F6.2, 1H)) 


WRITE OUTPUT TAPE 4,887,_,WCRK» ENERGT 


FORMAT(25 XSHWORK2E12.4, ICX7HENERGY=E12.4/) 
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20 


21 


22 


25 


24 
29 


30 


31 


32 


35 


STOP 


END 


SUBROLTINE STRESS (EPSI,EPSIA,SIGA,SIG) 


COMMON A,SIGY 
IF(EPST)20520,21 ; 
SIG=SIGA+(EPSI-EPSIA) 

GO TO 33 

IF (EPSI-EPSIA)22, 22,23 
SIG=SIGA+(EPSI-~EPSIA) 

GO TO 33 
SIGB=SIGA+(EPSI-EPSTIA) 
B=3.0/(7.0#(SIGY# #(A-1.0))) 


SIGI=EPSI 


SIGZ=SIGi—-(SIG1+B #(SIGI##A)-EPST)/( AeBe(SIGI#e(A-1.0))41.0) 


IFCABSF(SIG2#( 10. 0##(-8.0)) )-ABSF( SIG2-SIG1)) 29.30, 30 


SIGI=SI1G2 

GO TO 24 

IF (SIG2-S1GB) 3153 1532 
SIG=S1G2 

GO TO 33 

SIG=SIGB 

RETURN 

END 


END 
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APPENDIX V 


Maximum Force, Residual Deflection and Maximum Stress Data 

Tabulated on the following pages are the maximum force, 
residual deflection, and maximum stress data for tne three 
slenderness ratios and two eccentricity ratios for all values 
of ©< and © which were investigated for pulse loading. The 
values of maximum lateral deflection that were available are 
also included in the tabulation, These values were only availe 
able from the solutions to computer Program COL 1 used to obe 
tain graphical results. Thus, only a limited number of these 
values are included in the tabulation, 

The dimensionless time increments AT used for the numer- 
ical integration were the same for a given value of © in 
ail the solutions for both computer programs. These values 


are ag follows: 


6 Ace 
OvOr 0.10 
0.05 0.10 
0.10 0,08 
0.20 0,04 
O750 0,02 
1.00 O.01 


The values of AT are more than adequate to maintain the 
accuracy of the numerical scheme, The minimum number of 
increments for a natural period of first mode lateral vibra. 
tion of the column in the numerical process is 62. The values of 
WT were selected to display the complete history of tne 
dynamic response to the pulse Joading for the graphs shown 


in App. III using the maximum 900 points permitted in Sub- 


Vel 





routine GRAPH, 

The values of maximum stress tabulated on the following 
pages are maximum compressive stresses if there is no sign 
before the number, or maximum tensile stresses if there is 
a negative sigh before the number, In either case they are 
the maximum absolute values of stress that occurred in Spring 


l or Spring 2 of tne column model, 


Ve-2 








- ne — = 
. - = 
° 





| Maximum Residual 
Force Deflection 









[o.1o | 2.289 | 7.66x20 


; 
| 2.391 
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Maximum Force, Residual Deflection and Maximum Stress Data 


He2 © | L/p= 100 eR = 0.05 


Maximum Stress 1 | Stress e 

Deflection:| Maximum Maximum 

. (Ym/@) (Simax) (So max) 
0.953 | 2.1210" | .--- | 0.00282| 0.00074 
1.001 2624x107" is eo 0.00283 0.00078 
1121 | 2.5000" | =. | -o.buae7! 0.00100 
0.20 | 1.319 | 2.77x10* | --- - | 0.00289] -0.00129 
1.658 | 2.40x10% | ---.. | 0.00286] -0.00233 













® 


OL= 2,5 _ LQ- 100 e/ex 0.05 















Stress 1 | Stress 2 


fosos | 1.029 | 6uusno™ |---| o.oosas] 0.00055 

r7.65n0" |---| 0.00320] 0.00132 

[o.z0 | 2.470 | gumao® |---| 0.00327] -0,00n8e 
-2 

~@2. 






~ 


® 


Residual Maximum 
Deflection | Deflection 


(y--/e) .(Y¥m/@) 
















_=-- | 0.00325] 0.00275 
__---__|_ 0.00329] -0.00263 








; ® 
eo 
e A s 
« $ 4 CI i 


9.00x10 


/Va3- 
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Maximum Force, Residual Deflection and Maximum Stress Data 


os L/e= 100 eR = 0.05 


Maximum | Residual Maximum | Stress 1 | Stress 2. 

Force {Deflection | Deflection:| Maximum | Maximun |: 
jo.o1 | 0.954 | rstuxio |---| _0,00334| -0.00110 
0.05 | 1.069 | 1.28x10" |. --- | _0.00337|_-0.00130 
(0.20 | 1.602 | 1.980 | --- __|_0,00352| 0.00236 
0.50 | 2.237. | 2.20e107" |__--- __|_ 0.00356] _-0.00310 
2.00 | 2.633 | 2.05x107 | --- - | 0.00360|_-0.00313) 















oe 3.5 Le= 100 ef= 0.05 





















Residual 


Stress 1 | Stress 2 
Deflection | Deflection 
(y-/e) 


o.or | 0.955 | 1eieio* |---| 0.003491 _-0.00131 
0.05 | 1.120 | 2.0510" | --- "| 0.003541 ~0.00255 
Lowo | a.s1e | 2urac |---| 0,00360|_-0.00195 
0.20 | 1.72% | 3.2510 |---| _0.00370|_-0.00268 
pee 
a 










0.50 | 2ehu7 | 3483x107 
-" 





Veh 
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Maximum Force, Residual Deflection and Maximum Stress Data 















Ke i L/es 100. eR = 0.05 

| Maximum | Residual Maximum | Stress 1 | Stress 2 | 
Force {Deflection | Deflection | Maximum | Maximum 

0.05 | r1u7 | 2.8910" |---| 0.00366 | -0.00175] 

0.20 | 1.836 | 4.5610" |---| 0.00383 | -0.00261 








eel bee 70 ofa 0.05 

















Stress ] | Stress 2 


-o,or | 0.81 | 3.27e10"*|  --- __| 0.00315 | 000149 
0,05 | 0.840 | 3.1910*| ---__|0.00315 | 0.00150 
0.10 | 0.847 | 3.13x10"*| ---__| 0.00315 | 0.00253 


Residual . Maximum 
Deflection | Deflection 


(y-/e) ( ¥m/@) 













0.20 | 0.868 | 4.23x10*|  ---__| 0.00324 | 0.00263 
0.50 | 0.934 | 2.83x07*|  ---__| 0.00312 | 0.00190 
1.00 | 0.982 | 6.2mx0"*]  --- | 0.00272 


v=5 








all 


0.00224 











wares 


air 





[9-20 | 0,898 | 2-t0xa0""_ |---| 0.00381 | 0.00264 
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Maximum Force, Residual Deflection and Maximum Stress Data 


Ka 1.5 L/p= 70 eR = 0.05 


Maximum Residual Maximum Stress 1 | Stress 2 
Force Deflection } Deflection | Maximum Maximum 
(p/P) | (ye/e) ¢ |} (ime) | (Simax) | (Sama 
0.841 | 2.10x10"' | --- | 0.00381 | 0.00149 
0.854 2.09x107! Lovee ee 0.00381 0.00153 











0.20 | 0.987 | 2.100" |---| 0.00382 | 0.00184_ 
0.50] 1.206 | 204x107" |---| 0.00382 |-0.00260_ 






- ois 2 Le 70 | e/R= 0.05 




















Residual 
Deflection | Deflection 


( y-/e) 
=| 











Stress 1 | Stress 2 

( 
3.790 | __+--__| 0.00404 | 0.00249 | 
0.05 | 0.878 | 3.86x0' | __---.__| 0.00405 | 0.00158 
010 | 0.954 | 3.980" | = | 0.00106 | 0.00176 
0.20] 1.08 | 4.05m07 | ---. | 0.00x07 |-0.00228 
0.50 | 2.377 | 3.430" |---| 0.00403 |-0,00368 


0.00402 | 0.00368 | 





V~6 





oe 


[0.10 | 0.995 | 5.770" |---| 0.00422 | -0.00202_ 


APPENDIX V 


Maximum Force, Residual Deflection and. Maximum Stress Data 


Ka 2.5 L/e= 70 R= 0.05 


Maxinum Residual . Maximum Stress 1] | Stress 2 

Force Deflection } Deflection Maximum Maximum 

(P,/P. ) (y/o) (Ym/O) (Simax) (S. max) 

0.841 | 5.36x10"' | --- | 0.00418} 0.00149 
0.897 5652x107 eee eo 0.00419 | 0.00161 










0.50 | 1.489 | 5.03x10"" |---| 0.00422 | -0.00413 
/1.00 | 1.619 | 3.060! | --- - | 0.00422] 0.00393] 






- o&<s ] = 5m Le= 50 e/ex 0.05 













Stress 1] | Stress a 


) 
[=| 20007 | 0.00298 
Te] ex001a8 | 0.00240 

Tonto [0.6 | 2.7oao*| | o.conis} 0.0005) 
we 
or 
| 


Residual 
Deflection | Deflection 


(y-/e), 























[ 
0.00396 
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Maximum Force, Residual Deflection and Maximum Stress Data 


HK 1.5 7 ae (R= 0.0 








; Maximum Residual ' Maximum Stress 1 | Stress 2 
| Force {Deflection Deflection | Maximum | Maximum 
ee a — a a = 





O. | 0.01 O. 0.585 5. ls O. 0.00446 0. | 0.00234 | 


[0.05 | 0.612 | 5.16x10"" |---| 0.00446 | 0.00252] 
[0.10 | 0.654 | ‘5.13m0-' |---| 0.00446 | 0.00275) 
ean [ oro a“ ed ee mer eco a 
(0.50 | 0.822 | 403x107! |---| 0.00445 | 0.00383) 
Cee eC 











ot = 2 | a oe 50 bai 0.05 


| Sait cia Residual Maximum | Stress 1| Stress 2 
Force | Deflection | Deflection | Maximum | . Maximum 
ies 
-o.or | 0.585 | 7.2710 |---| 0.00461] 0.00235 
rov0s | 0.626 | rasa" |---| o,0046e| 0.00057 
[eae [exes | ranact | [ eoot scone] 
on20 | 0.756 | 726ao" | === | 0.00nee| 0.00938 
Fo.s0 | 0.473 | aga" | --- 0.00461 | -0.00%54 
[avo [ons | zone" [ --- | 0.00136) 0.00139 
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Maximum Force, Residual Deflection and Maximum Stress Data 


XL « 205 | Lee 50 OR «= 0.05 


— Residual : Maximum | Stress 1 | Stress e 

Force {Deflection | Deflection | Maximum } Maximum 

Tovar | 0.586 | oaumo™ |---| o.o0ure | 0.00036 
ro.05 | 0.645 | ozone" |---| 0.00073 | 0.00269 
[0.10 | 0.700 | 9.4ox0" |---| 0.00473 | 0.00300 | 
0.20 | 0.794 | 8.980" |---| 0.00474 | -0.00013. 
10.50 | 0.909 | 6.17x10"'\ |---| 0.00475 | -0.00475. 
[2,00 | 0.948 | 2zaxro |---| 0.00471 | 0.00456. 








o= 2 Le 100 e= 0.01 





















Residual Maximum Stress 1 | Stress 2 | 


Deflection | Deflection 
(y-/e) (ym/@) 


0.01 | 1.003 | 2.s6xo* |---| 0.00287 | _0,00100 
0.05 | 1.192 | 3.03x0*| ---_|_ 0.00292 | 0.00118 
0.20 |.1.616 | 3.30x107*| 1.81 __| 0.00295 | -0.00180 
0.50 | 1.911 | 3.07x07*| 1.93 | 0.00294 | -0.00231 | 
1.00 | 1.98¢,|- 2.951073] 0.69 | 0.00250 | 0.00235 
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Loo | 2.578 | 2.0710 |---| 0.00354 | -0,00224 
[0.50 | 2.535 | z.1axi0% |---| 0.00359 | -0.00326 
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Maximum Force, Residual Deflection and Maximum Stress Data 


See L/p= 100 eR= 0.01 


Maximum Residual Maximum Stress 1 | Stress 2 
Force Deflection } Deflection faximum Maximum ‘|: 

(Pe. /%, ) (y-/e). —— —- 2 
0.01 Ol 1.017 6. 93x10 | mee @. | 0.00317 | OY | 0.00200 | 
1.258 | @.eox0% |---| 0.00325 | -0.00139 











0.20 | 2.813 | 127x107! | ---__|_0.00334 | -0,00250_ 
0.50 | 2.264 | retaxio~ |---| 0.00335 | 0.00286 
j1.00 | 2.423 | 2.zoa0t] --- | 0.00306 | 0.00287] 








o< = nl LA= 100 e/e= QO.O1 










Residual Stress 1 | Stress 2 


Deflection | Deflection 


(y-/e) 


wm & 
oor | 2,032 | 1.3010" |---| 0.00338 | -o,00115 
0.05 |" 1.320 | r.72eot |---| 0.00348 | -0.00173, 










0.20 | 1.971 | 247x107 |---| 0,00360 | -0,00274, 






1.00 | 2.7m | 7.27e0t| --- | 0.00346 | 0.00327 
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Maximum Force, Residual Deflection and Maximum Stress Data 


? 


Ho 3.5 Ifez= 100 . = *R» 0.01 


a Maximum | Residual Maximum | Stress 1 | Stress 2 
a Force {Deflection }| Deflection | Maximum | Maximum 
—{o,or | 1.048 | 2,00x0" |---| 0.00353 |-0.00135_ 
[0.05 | 1.380 | 2.6907 |---| 0.00363 | -0.00199_ 
| 3.26107 |---| 0.00370 |-0.00255_ 
0.20 | 2.108 | 3.86x1o" |---| 0.00377 |-0.00292 
0.50 | 2.739 | 3.o0x07 | --- _. | 0.00376 {-0.00358 | 
1.00 | 2.961 | 1.35m0% |---| 0-00372 | 0.00354 | 











oe 7 L= 100 ep= 0.01 



























Residual -_Maximum Stress 1 | Stress 2 
Deflection | Deflection 
(y-/e) (ym/@) 


ae ie | ene he eee 
oe | is | worl = (een ee 















0.20 | 2.230 | 5.u9xio | 2.97 | 0.00390 | -0.00308 


Vell. - 


=| 








a 





LE EEE 


[[oxos | o.9ar | a-eomo™] -—- | 0.00223 | 0.00001 


ar 
| =~ {o. 
Loto | 1.072 | 2.zzaot |---| 0.00383 _ 
a 2 Ee 
_ =~ |o. 
_ = fo. 


Maximum Force, Residual Deflection and Maximum Stress Data 


Heo. 1 | L/e= 70 eRa 0.01 


Maximum Residual Maximum Stress 1 | Stress 2 
Force Deflection Deflection:| Maximum Maximum 
(P/P. ), = /2) (Ym/Q) —- oo 






















o<= 1.5 


















~m e 
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Maximum Force, Residual Deflection and Maximum Stress Data 


Xn 2 L/e= 70 e-= 0.01 


Maximum Residual Maximum Stress l] | Stress 2 
Force Deflection } Deflection Maximum Maximum 
(227. ) (y-/e) ; (Y¥m/Q) (Simax) (Somax) 
0.930 | 4.oimio’ | --- | 0.00406 | 0.00200 
7 1as1 | y.azac |---| 0.00408 | 0.00250 
| 1.497 | 3.04x10"' | --- - | 0.00408 | 0.00359 
}1.00 | 1.542 5.20x107* oe 0.00388 0.00380 | 







= 2.5 Le= 70 eR= 0.01 













Residual Maximum Stress 1 |] Stress 2 


Deflection | Deflection 
(y-/e) (ym/@) 
















ae 





. 
® 
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Maximum Foree, Residual Deflection and Maximum Stress Data 


cs | L/e = 50 | eR = 


Residual | Maximum | Stress 1 | Stress 2 
Deflection } Deflection Maximum Maximum 
(¥r/Q) (Ym/Q) (Simax? | (Sama) 











. - “& 
Gee? \gones3 | zegazio"’ | --- 0.00419 | 0.00284 
mer oie ny 


0.10 | 0.689 0.67 | 0.00421 | 0.00313 
0.20 | 03730 265 x00-' 0.75 0.00421 | 0.00337 


0.50 | 0.775 | 1.15x107' | 0.48 0.00403 | 0.00367 
1500 | Ouge2 | iepexio—* 0.10 





o< s LS Les 50 es Ose 
| Max4dnum Residual Maxinun Stress 41 | Stress 2 
Force Deflection | Deflection 
(P/F, ) (y¥r/@) (¥m/Q) 


0.01 | 0.636 | 5.30x107' | we= | 0.00447 | 0.00286 | 
| 0,05 | 0.687 | 5.22x107 pe ee 
| 0,10 | 0.735 | 5.02x107 eZ a 
reso [ease [pear | oe fem [om 


Veld 





a as ee 
‘ « 
. 
. - e 
e 


Maximum Residual Stress 1 | Stress 2 
Force Deflection Maximum |° Maximum 
| (P/% ) (y-/e) . Corea (Soe aay) 


[aco [over | 3.66n0 | one 
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Maximum Force, Residual Deflection and Maximum Stress Data 


Ks 1 L/e= 50 | = (eRa 0.01 










ae 
Deflection 
(Ym/@Q) 
recor | 0.63 | ammo" | | o.oo1s | 0,0028s 
0.05 | 0.658 | a.rmao™ |---| 0,00u20 | 0.00097 | 






Se 1.5 | L= 50 efea 0.01 














Stress 1 | Stress 2 


Residual 
Deflection | Deflection 


(y-/e) 
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Maximum Force, Residual Deflection and Maximum Stress Data 


A= 2 | L/e= 50 ek-= 0.01 


Maximum | Residual Maximum Stress l Stress 2 
Force {Deflection } Deflection ‘faximum | Maximum | . 
Tawar" |---| 0.00463 | 0.00289 | 
To.0s | 0.709 | rorwao™ |---| 0.00463 | 0.00323 _ 
Tonio] 0.766 | rarao" |---| 0.00464 | 0.00354 
o.20| 0.832 | 6.27mom |---| 0.00463 | 0.00991_ 
5,50 | 0.907 | aeauino~ | _--- _| 0.00463 | 0.00435_ 
boos | ucbuxio® |---| 0.00049 | 0.0006 | 





































o< = 2.5 | L/Q= 50 ees 0,01 
Residual Maximum 
Deflection | Deflection 


7 Maximum Stress 1 | Stress 2 

Force Maximum | - Maximum 
8 (P,/Pe) |  (y¥r/e) (¥m/@) (Semax) 
5.25x10-2| 0422 | 0.00465 | 0.00462] 
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APPZNDIX VI 

Evaluation of the Approximations and of the Accuracy of 
the Numerical Scheme 

From a consideration of the maximum lateral deflections 
encountered in the solutions the approximations made in de- 
veloping the equations for the column model may be evaluated, 
The maximum value of © for all solutions was less than 0.08 
radian, Based on a maximum © = 0.08 radian the errors are 
as follows: 


(a) 0.01 per cent error in (L/2-e tan €}) as a result 
of neglecting e tan® 


(b) 0.2 per cent error in [-t-sin@-e(1-cos-@)| by con- 
sidering sim@=€© and neglecting e(l-cos©@) 


(c) 0.25 per cent error in tam@and 0.11 per cent 
error in sin® as a result of considering 
tan @ =sine=ecC 


(ad) less than 0.06 per cent error in (l-cos©@) as a 
result of considering l-cos@=4-4* 


From this tne approximations which have been made in deriv- 
ing the equations for the column model are reasonable and 
produce insignificant error in the results, 

As an additional check for cumulative error a compar- 
ison of the ram work done on the column model witn the net 
energy input to each of the springs (a combination of strain 
energy of plastic deformation and the existing elastic strain 
energy) is made in all computer solutions for the results 
shown in Figs, 10 thru 14, These agree within 1 per cent 


thus verifying the compatibility of the equations and check-— 





ing the accuracy of the numerical scheme for the time incre- 
ments chosen, Work-energy considerations are not included 
in the computer program for obtaining grapnical results but 
the dimensionless time increment AT used for the numerical 
integration in both programs for a given value of © was the 


same, 
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